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Abstract 

The open light gluino window allows non-trivial higher twist gluino contribu- 
tions to the proton wave function. Using a two-component model originally 
developed for charm hadroproduction, higher twist intrinsic gluino contribu- 
tions to final state i?-hadron formation are shown to enhance leading twist 
production in the forward xp region. We calculate i?-hadron production at 
Piab = 800 GeV in pp, pBe, and pCu interactions with light gluino masses of 
1.2, 1.5, 3.5, and 5.0 GeV. 

PACS Numbers: 12.60.Jv, 13.85.-t, 14.80.Ly 
INTRODUCTION 

The gluino is the supersymmetric partner of the gluon. It is an electromagnetically 
neutral, adjoint fermion with the same color structure as its boson counterpart. As yet, no 
clear experimental evidence of supersymmetric particles has been found. The most likely 
reason for this is the large expected mass of the supersymmetric particles (Asusy ~ 1 
TeV). However, an intriguing scenario exists whereby the gluino is not only the lightest 
supersymmetric particle but also very light compared to the SUSY scale, <^ 100 GeV. 
This possibility arises naturally in a number of quite attractive models characterized by 
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special boundary conditions at the grand unification scale [||,|| and in certain models of 
gauge-mediated supersymmetry breaking 

Light gluinos are predicted to form relatively light bound states of quarks or gluons and 
gluinos called i?-hadrons [Q. The lightest predicted i?-hadrons include mesinos {qqg), two 
barinos, R^{uudg) and S^{udsg), gluinoballs (gg), and the glueballino or (gg)- The 
properties of i?-hadrons including their mass, decay modes, and lifetimes depend strongly 
on the mass of the gluino. 

There have been many theoretical and experimental attempts to find evidence for and/or 
exclude the light gluino scenario. Searches for i?-hadrons produced in fixed target experi- 
ments have been performed for a number of the predicted i?-hadron decay channels 
Effects of a light gluino on QCD observables have been analyzed [^. Stable particle searches, 
T decays, beam dump experiments etc. all have potential sensitivity to the presence of a 
light gluino or the i?-hadrons. A brief summary of the various possible resulting constraints 
on a light gluino is given in Ref. |T^. In addition, Ref. [rT| claims that rrig > 2.5 — 3 GeV is 
excluded on the basis of their analysis of OPAL data. Although these various analyses are, 
in combination, potentially sensitive to most regions of light gluino mass, all rely on model- 
dependent inputs. As a result, we believe that at present it is impossible to definitively 
exclude any gluino mass below 4 — 5 GeV. Thus, it is of great interest to find additional 
approaches for discovering and/or constraining light gluinos and the i?-hadrons. 

In this paper, we will explore the possibility of detecting i?-hadrons at large xp in pp 
and pA fixed-target interactions. Our calculations will be restricted to the rrig ~ 1.2 — 5 
GeV region where we can be confident that the semi-perturbative techniques that we employ 
are reliable. This region is of particular phenomenological interest because of the analogy 
that can be drawn between heavy quark and light gluino production. Indeed, if the gluino 
and heavy quark masses are comparable, one might anticipate observation of hard gluino 
production analogous to that already observed in high-Xi;' charm hadroproduction ||12[. The 
leading-twist pQCD predictions for charm production in pp and pA collisions fail to account 
for many features of the high-Xi? data. These include unexpectedly large production rates 
and anomalies such as flavor correlations between the produced hadrons and the valence 
spectators, manifested as leading charm and a strong / asymmetry in Tr~A interactions 



||12|| , double J/\t' production at large xp |]13[) ^^'^ Feynman scaling of J/^ production in pA 
interactions all of which suggest a breakdown of factorization |]T6[ at large xp- The 

anomalies and cross section enhancement may be partly explained by higher twist terms in 
the operator product expansion (OPE) on the light cone associated with the dynamics of 
the QCD bound state. Analgous terms should be present for light gluinos. 

The intrinsic charm model (IC) |]T7|>II11 approximates non-perturbative higher twist Fock- 
state contributions of heavy quarks in hadronic wave functions. The phenomenological 
predictions of IC directly address the above puzzles in charm hadroproduction |ll9|-p4|. For 
example, IC provides a coalescence mechanism whereby final state hadrons can share valence 
quarks with the projectile, naturally producing leading particles. 

In analogy with leading charm, we study i?-hadron distributions using "intrinsic gluinos" 
(IG) in regions of phase space where the gluino mass and momentum fractions conspire so 
that higher twist effects cannot be ignored. In this paper, we calculate enhancements over 
the leading twist i?-hadrons xp distributions with gluino masses rrig = 1.2, 1.5, 3.5, and 5.0 
GeV. Both pp and pA interactions at pi^h = 800 GeV are considered. 
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PQCD LIGHT GLUINO HADROPRODUCTION 



In pQCD, gluinos are produced in pairs by gg fusion and qq annihilation, gg, qq — > gg, 
as well as quark-gluon scattering to squark and gluino, qg qg. Precision Z-pole data has 
constrained the squark mass to be greater than 100 GeV, quite large compared to the light 
gluino masses considered here. Therefore, we expect that the qg contribution with the virtual 
squark in the t-channel will be small compared to the other contributions, particularly at 
fixed-target energies. 

The leading twist inclusive i?-hadron xp distribution at leading order is obtained from 
the gluino xp distribution {xp = (2mr/ys) sinh?/) which has the factorized form in pQCD 



da ^/s 
dx^ ~ "2~ 

^ ^,3 



dz2,dy2d^PT 



El 



Z3 



VT dt 



(1) 



Here a and b are the initial partons from projectile and target hadrons A and B, 1 and 2 are 
the produced gluinos, and 3 is the final-state i?-hadron. The sum over i and j extends over all 
partonic gluino production subprocesses. A K factor of 2.5 is included to account for NLO 
corrections. Since the K factor is approximately constant with xp for charm production 
except as Xi? — 1, we assume that the K factor for gluino production is also independent 

of Xp. 

The fragmentation functions, DH/g{z) with z = Xn/xg, describe the collinear fragmen- 
tation of final state i?-hadrons from the produced gluinos. For simplicity, a delta function 
was used for hadronization, Dn/giz) = 6{z — 1). This assumption results in the hardest xp 
distribution at leading twist since the i?-hadron carries all of the gluino's momentum. Other 
fragmentation functions would soften these distributions. Note that for any fragmentation 
function to factorize, it must be independent of the initial state {i.e. it only depends on 
z^ and not Xa). Thus, regardless of the fragmentation function used, all i?-hadrons will be 
decoupled from the initial state to leading twist. 

The partonic cross sections for gluino production in Eq. ([l|) are ||25[ 
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(s — m? + m?)(m? — i) — m^s 1 s{u + m?) + 2(m? — m?)(m? — m) 
s(m| - 1) 18 s{u- mf) 

1 (m? - i) (2m + m? + t) 1 (m? - t) (s + 2t - 2m?) 
^ 2 2(£-m?)(n-m?) ^ 2 2(£ - m?)(M - m|) 
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We calculate leading twist pQCD gluino distributions for 800 GeV pp interactions. Fig- 
ure shows the gluino distributions using the MRS D-' parton distributions in the proton 
||26| with rrig = 1.2,1.5,3.5, and 5.0 GeV and rriq = 100 GeV. The characteristic falloff at 
large is similar to heavy quark production. Choosing a larger squark mass would only 
marginally decrease the total cross section because the gq channel is suppressed by the large 
squark mass. The gluino production cross section is a strong function of mass. The cross 
section is largest for rrig = 1.2 GeV and decreases by a factor of 3 for rrig = 1.5 GeV. There 
is then a drop of 250 to the rrig = 3.5 GeV gluino cross section and another factor of 20 
between the 3.5 and 5 GeV cross sections. Additionally, the falloff of the cross section with 
xp becomes steeper as rrig is increased. 

Charm hadropro duct ion phenomenology has taught us that higher twist contributions 
can become comparable to leading twist in certain parts of phase space, introducing cor- 
relations between the initial and final states. These effects will be addressed in the next 
section. 



INTRINSIC CONTRIBUTION TO HIGHER TWIST 

In deep inelastic scattering, higher twist terms in the OPE are suppressed by a factor of 
l/Q^"". These terms are essentially irrelevant when is large. Analogously, in hadropro- 
duction, a similar suppression of 1/M^ typically renders higher- twist effects unimportant 
except in regions where pQCD is seemingly inapplicable {i.e. where is small). However, 
it has been shown that in the simultaneous oo and x — 1 limit with M^(l — x) 

fixed, a new hard scale emerges where higher twist contributions to the cross section be- 
come comparable to leading twist P7|-PU|. In the case of heavy quark production, this new 
scale can be associated with either the resolution of the transverse size of the intrinsic heavy 
quark pair or with the transverse resolution of any "pre-coalesced" hadrons inside the parent 
hadron. The heavy quark fluctuations can carry a large fraction of the projectile's forward 
momentum since the constituents of the bound state move with the same velocity. The Fock 
state may be broken up by an interaction with soft gluons in the target, producing a leading 
hadron containing a heavy parton. 

The bound state wave function for a state containing higher-twist contributions can be 
obtained from the Bethe-Saltpeter formalism evaluated at equal "time" on the light cone 



17,30 



{M' - J^ti'^MxukT.) = j\dy] J ^^K{x,,krr,y,Jri,M'MyM (5) 
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where M is the mass of the projectile hadron. The transverse mass of an individual parton 
is defined by rhi^ = k^, + mf, where fc^i is the transverse momentum of the i^^ parton in 
the n-particle Fock state, \qi, g^, The momentum fraction of the i^^ parton in the 
Fock state is Xj, [dy] = n"^^(iyj5(l — S^^j^y,) is a longitudinal momentum conserving metric 
and [d'^lr] = n"^;^(i^%5^(S"^;^/jv). The interaction kernel is K. 

The simplest way to create final state hadron distributions from a specific Fock state 
wave function is now described. The vertex function on the right hand side of Eq. is 
assumed to be slowly varying with momentum. The operator on the left hand side of the 
equation is then evaluated at the average transverse momentum of each parton, (k^.), with 
the constraint ^2^Ti = 0. With these assumptions, the transverse mass of each parton 
is fixed and the vertex function becomes constant. The probability distribution is then 
proportional to the square of the wave function which is now inversely proportional to the 
off-shell parameter A = — S"^]^ ('^i ) / where {fhf) is the average transverse mass squared 
of the i^^ parton. After longitudinal momentum conservation is specified by 6{1 — 
the probability distribution becomes 



f f 2 dx ^ 



Nn6{l - Sr^ixOA-^ (6) 



where Nn is the normalization constant for an n-particle distribution. The probability 
distributions as a function of x for any final state hadron can be generated by integrating 
Eq. (H) including final state coalescence constraints. 

The characteristic shape of the longitudinal momentum distribution of the final state 
hadron can now be obtained up to an overall normalization constant. The important feature 
of this model is that final state particles are not "produced" in a collision, as such, but are 
rather "intrinsic" to the projectile's Fock state and are liberated after a soft interaction 
with the target. This intrinsic source of final state particles acts as a perturbation to the 
dominant parton fusion mechanism. However, unlike parton fusion, it incorporates flavor 
correlations between the initial and final states. This mechanism will dominate the total 
cross section in the limit xp —>■ ^ since xp ^ x when the final-state hadron evolves directly 
from the projectile wave function. 

In this paper, we assume that the model developed for heavy quark hadroproduction at 
higher twist can be applied to gluino production in the proton wave function. Final-state 
i?-hadron production from IG states is described in the remainder of this section along with 
its relationship to IC production. The characteristic shapes of the intrinsic distributions in 
the proton were generated for the gluino alone and for the R^{uudg) and S^{udsg), and 
the R^{gg). In all cases, the "minimal Fock state" was used to generate the final state 
coalescence. This emphasizes the most leading final states. 

The gluino can fragment into a i?-hadron, just as in pQCD production. In this uncorre- 
lated case the hadron xp distribution is 



pkF ^ D I i 7\ 

-j^ = NkJ U^.^,dxjdz5{l -E^») J' ^ (^g - zXg)A-'' , (7) 



where k indicates the order of the Fock state containing the intrinsic gluinos {i.e. the x^'s 
are included among the Xi). Gluinos are produced in pairs because other supersymmetric 



5 



vertices involving squarks and photinos are highly suppressed due to their much greater 
masses. The minimal proton Fock state with a gluino pair then has five particles, \uudgg). 
Fragmentation of other, higher, Fock states will have a smaller production probability and 
produce gluinos with lower average momentum. For consistency with the coalescence pro- 
duction described below, we include fragmentation of six and seven particle Fock states with 
R? and production respectively. 

i?-hadron production by coalescence is specific to each hadron. The intrinsic gluino Fock 
states are fragile and can easily collapse into a new hadronic state through a soft interaction 
with the target, as is the case for IC states. The coalescence function is assumed to be 
a delta function. The momentum fraction of the of the final state hadron is the sum of 
the momentum fractions of the of the i?-hadron valence constituents from the proton wave 
function. The three i?-hadrons we consider are all calculated from only the minimal Fock 
state required for their production by coalescence. Thus, only the most leading configuration 
is used. As in the fragmentation case in Eq. (|^, including higher Fock components does not 
significantly increase the total rate because the other Fock state probabilities are smaller 
and also does not enhance the yield at large xp because the average xp oi coalescence is 
reduced relative to that from the minimal Fock state. 

The five-particle Fock state \uudgg) produces the most leading i?-hadron, the i?^, be- 
cause the is generated from four of the five constituents of the Fock state. 

= iVsP^ j Y{]^^dxj5{l - Xi)5{xR+ -Xu-Xu-x^- x^) A-^. (8) 



dxji^ 



i=l 



Here, is a factor incorporating the coalescence probability given the five-constituent Fock 
state. Note that in this case, the i?"*" xp distribution is proportional to the gluino distribution 
in Eq. (|^), obtained by setting DH/g{z) = 6{1 — z), with k = 5 evaluated at 1 — xp- 

The is generated from a six-particle Fock state, \uudggg). Unlike the gluinos, single 
gluons can be included in the higher-twist Fock state since one gluon can couple to two 
quarks in the Fock state. The six-particle state is the most leading state for R^ production. 
The coalescence of R^ hadrons is described by 

dXRO J .^^ 

The last i?-hadron we consider is the which, since it contains an s quark, must 
be produced from a seven-particle Fock state, \uudssgg). The will have a harder xp 
distribution than the even though the average momentum fraction of each constituent 
in the seven-particle state is smaller than those of the six-particle state. This harder xp 
distribution is due to the greater number of constituents, four, rather than the two 
constituents. In this case, 

*9 — ^^^p7 



dxso 



N-jPl j Il]^^dx,5{l - Xi)5{xs, -Xu-Xd-Xs- Xg)A-^. (10) 

i=l 



In what follows, the coalescence probabilities P^'^'^ appearing in Eqs. (|), , (p!0D are 
taken to be unity. That is, it is assumed that the gluinos will always coalesce. 
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Figure |^ shows (using arbitrary normalization) the characteristic x dependence of the 
probabihty distributions in Eqs. (|7|)- (|T0|) with rrig = 1.2 GeV. The single gluino distribu- 
tion is calculated using k = 5 and Dn/giz) = 6{1 — z) in Eq. (J^). _R-hadrons produced 
by uncorrelated fragmentation have the softest xp distributions, {xg) = 0.24 when k = 5. 
Contributions from progressively higher single gluino Fock states have smaller relative prob- 
abilities, as we discuss below, and a decreased {xg), which would eventually build up a gluino 
sea in the proton. The distributions from coalescence are all forward of the single gluino 
distribution. As expected, since the takes all three of the proton valence quarks, it is 
the most leading i?-hadron with = 0.76. The distributions for the other final state 

particles, the 5*0 and the Rq, are softer with (X50) = 0.56 and {xj^o) = 0.35 respectively. 

We have shown the results with the lowest gluino mass we consider. Increasing the mass 
increases the average xp of the gluino distribution of uncorrelated fragmentation, Eq. (^, 
but leaves the average xp of the R hadrons unchanged in the mass range we consider. 

The intrinsic gluino production cross section for i?-hadrons, from an n-particle Fock state 
is written by analogy with the IC cross section 

'^Upp) = GcP:gat{mgg)a^^^ , (11) 

where Gc is a color factor. The inelastic pp cross section is ~ 35 mb at 800 GeV. The ratio 
/Arh?g sets the scale at which the higher and leading twist contributions are comparable. 
We use /i^ ~ 0.2 GeV^, consistent with attributing the diffractive fraction of the total J /if) 
production cross section to IC |T9|J20| , |2^ . There is a factor of because the intrinsic state 



couples to two of the projectile valence quarks. The higher-twist contribution then contains 
two more powers of than the leading-twist contribution. This factor is included in the 
cross section rather than in the probability distributions as done previously |^,24| to more 



explicitly show the effect of this dependence on the cross section when the mass of the 
intrinsic state is changed. 



Since the intrinsic charm cross section is 
the two cross sections are related by 



{PP) = P>tim.Ml7^-. ' (12) 



cr?^iPP) _GcPrsrhlat{mgg) 



(^Upp) P^c rhlaiirricc) 



(13) 



The relative color factor between intrinsic gluinos and intrinsic charm, represented by Gc, 
may enhance the IG contribution over that of IC because of the color octet nature of the 
gluino. However, in this work, to isolate mass effects, we assume the color factors for IG 
are the same as IC, setting Gc = 1- Changing Gc would effectively scale the cross section 
ratio in Eq. (|13D by a constant factor. The overall effect of changing Gc is small relative to 
the leading-twist cross section unless Gc is very large. The intrinsic charm mass is used as 
the scale from which to approximately evolve the intrinsic gluino cross section as previously 
done for intrinsic beauty Note that when Gc = 1, if rhg = ttIc, the IG and IC cross 
sections are the same. The IG cross sections are normalized by scaling Pig in proportion to 
Pic, as described below. 
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A limit of Pj^ = 0.31% was placed on the intrinsic charm probability in the five- 
particle state \uudcc) by charm structure function data p2|- p^] . The higher Fock state 
probabilities were obtained from an estimate of double J/\l/ production [0, resulting in 
^icc ~ 4.4%P^^ [^. Mass scaling was used to obtain the mixed intrinsic charm probabili- 
ties, Pjg^ ~ {rhc/rhqYPj^^ |T^. To obtain the n-particle gluino Fock state probabilities, Pj^, 



we assume that the same relationships hold for the gluino states. The five-particle gluino 
state then scales as 

" 1 

Tf) 

4 = -:%Pl ■ (14) 
^ rrig 

Assuming P^^^ = 4.4%Pj^, the seven-particle Fock state probabilities are 

Thus, if rhg = rhc, Pf^ = Pfg and Pj^^ = Pj^g. For simplicity, the probability for the 
mixed gluon-gluino proton six-particle Fock state was set equal to the seven-particle mixed 
probability with rhg = rhg. The effective transverse masses used were rhg = rhg = 0.45 GeV, 
Trig = 0.71 GeV, and rhc = 1-8 GeV. The transverse mass of the gluino, rhg, is fixed to the 
values of rrig used in the leading twist calculation. 



COMPOSITE MODEL PREDICTIONS 

In this section, we calculate the total xp distribution of final-state P-hadrons including 
both leading- and higher- twist contributions. The model predictions for P"*", P" and 
production on proton and nuclear targets are then given at 800 GeV. 

The final state da/dxp distribution is the sum of the leading twist pQCD distribution 
and the higher twist intrinsic contributions. Since many experiments use a nuclear target, 
the characteristic A dependence of each contribution is included, 

A—^ + A^—^ . (16) 



dxp dxp dxi 

The first term is the leading twist term whereas the second term is the higher twist IG 
contribution. Leading twist necessarily involves single parton interactions between the target 
and the projectile and thus cannot account for collective nuclear effects. Thus, the leading 
twist cross section scales linearly with the number of nucleons in the target modulo nuclear 
shadowing effects. The nuclear dependence of J/ifj production in pA interactions shows 
that if the nuclear dependence is parameterized by A"", a— >2/3asxi;'— i>l [0,|1^]. The 
emergence of this surface effect at large xp is consistent with spectators in the projectile 
coupling to soft gluons from the front face of the target rather than the volume. The 
NA3 collaboration extracted the A dependence of J /ip production at large xp and obtained 
j3 = 0.71 in Eq. (|TB|) \T^\. We use the same value of (3 for charm production since the 



available data on the charm A dependence [|35| leads us to expect a similar A dependence 
for charm and J/il^ production at large xp. 
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The intrinsic gluino contribution to i?-hadron production includes contributions from 
both hadronization of single gluinos by uncorrelated fragmentation, Eq. and coalescence 
into final-state i?-hadrons, described in Eqs. (^-(|1^). That is, 

= ^1^ + ^2-7^ 17 

dxp dxp dxp 

where P"^^ and P^^ are the IG contributions from fragmentation and coalescence respec- 
tively. The parameters and .^2 allow adjustment of the relative gluino fragmentation and 
coalescence contributions. We used single gluino fragmentation from the same Fock state 
as the coalesced hadron. That is, for i?"*". A; = 5 in Eq. (|^, while k = 6 for and k = 7 
for S^. We fix G = ^2 = 0.5 for simplicity. For a more realistic accounting of all possi- 



ble contributions to Eq. ([T7|) for charm production, see Ref. for relative charm hadron 
production probabilities in the proton. The respective fragmentation and coalescence prob- 
ability distributions in Eq. ( p!?] ) are converted to cross sections using Eq. (|lT]) and added to 
the leading twist cross section as in Eq. ([T6|) . 

We calculate _R-hadron production at 800 GeV in pp, pBe, and pCn interactions with 
TUg = rhg = 1.2, 1.5, 3.5, and 5.0 GeV. Delta function fragmentation was used for single in- 
trinsic gluino production by uncorrelated fragmentation and for leading twist hadronization. 
That is, we take Dh/q^z) = 6{1 — z) in Eqs. (|I]) and (0). 

Figure 1^ shows the normalized i?-hadron xp distributions calculated according to Eq. ([TlD 
in pp interactions with = 1.2 GeV. The difference in the yields as — >■ is due to the 
difference in probability for the five, six, and seven particle Fock states. The and 
cross sections are similar at low xp because we have assumed Pf„^ = Pj„^, as described in 
the previous section. However, the shapes are different at low xp because the probability 
distribution for uncorrelated fragmentation has a smaller average {xp) when = 7 in Eq. (^. 
The has the largest cross section of the three hadrons. Its distribution is symmetric 
around xp = 0.5 because the fragmentation yield and the R'^ yield from coalescence are 
symmetric in the five particle Fock state. The yield increases near xp ~ 0.25 due to 
the forward peak of the coalescence distribution seen in Fig. |[ The yield at low xp is 
relatively reduced because the fragmentation calculation with = 7 is narrower so that the 
two peaks are effectively separated in Fig. ^ Since the fragmentation peak for k = 6 and 
the coalescence distribution lie close together, they blend into a broad peak for the R^ 
Xp distribution. 

Figures ^, ^, and |^ show the predicted i?"*", S^, and RP xp distributions per nucleon in 
pp, pBe, and pCu interactions at 800 GeV calculated according to Eq. (|T6|). Each figure 
includes all four gluino masses. As Xi? — the xp distributions of all targets are equal for a 
given nig. This indicates the dominance of leading twist production at low xp, independent 
of the final state. As xp ^ 1 the higher twist terms begin to contribute. These higher 
twist effects are suppressed in nuclear targets because of their slower relative growth as a 
function of A compared to the leading twist A dependence. Although larger mass gluinos 
are more difficult to create, the relative contribution to the total cross section from higher- 
twist production in Eq. (pIBf ) increases with gluino mass because of the slower decrease of 
the intrinsic gluino contribution relative to the mass suppression of the leading twist cross 
section. The greater mass suppression of the leading twist cross section also influences the 
value of Xp where the higher twist contribution begins to appear. Increasing the gluino mass 
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leads to intrinsic gluino effects appearing at lower xp- This effect is seen in Figs. When 
rrig = 1.2, IG effects become obvious near xp ~ 0.5 while IG contributions begin to appear 
for xp ~ 0.2 in production when rrig = 5.0 GeV. 

Dramatic leading effects are predicted for the which, as pointed out above, shares 
three valence quarks with the proton in a minimal five-particle Fock state configuration. 
This characteristic "hardening" of the xp distribution for xp > 0.6 should be clear in a 
successful i?"*" search. However, the leading effects are also present for the other particles. 
The is the next hardest distribution, sharing two valence quarks with the proton while 
the tends to be the softest, since no projectile valence quarks are shared. 

For a clearer comparison of the leading effects predicted for each final state i?-hadron. 
Figs. show the and distribut ions together in pp interactions with rrig = 

1.2,1.5,3.5, and 5.0 GeV respectively. The leading twist gluino distribution is also shown 
for comparison. In each case, the intrinsic contribution begins to emerge from the leading 
twist calculation between xp ~ 0.2 and xp ~ 0.4. In Fig. ^ with rrig = 1.5 GeV, the 
predicted i?"*" enhancement at xp ~ 0.8 is about 700 times larger than the leading twist 
prediction. At the same value of rrig and xp, the contribution is about 40 times greater 
while the R^ is just under 6 times greater. When the gluino mass is increased to rrig = 5.0 
GeV, shown in Fig. ^ the i?" dominates i?-hadron yields for xp < 0.6. This is a consequence 
of the increased (xp) for single gluino fragmentation at the larger mass. Although the cross 
sections are small at rrig = 5.0 GeV since the gluino mass is comparable to the bottom mass, 
the predicted enhancements over the leading-twist baseline are quite large: 2.5 x 10^ for the 
-R"*", 1.6 X 10^ for the S^, and 281 for the The enhancements are in fact larger than those 
with smaller gluino masses due to the greater mass suppression of the leading twist cross 
section. 

CONCLUSIONS 

The light gluino window opens the possibility of non-trivial higher twist gluino contri- 
butions to the proton wave function. In analogy to charm hadropro duct ion, intrinsic gluino 
Fock components contribute to final state i?-hadron formation, enhancing gluino production 
over leading twist parton fusion in the forward xp region. 

In this work, we have studied a "maximally leading" scenario for final state i?-hadrons 
in pp and pA interactions at 800 GeV. Our model predicts that the contributions of higher- 
twist intrinsic states lead to strong flavor correlations between initial and final states for 
Xp > 0.6. The large intrinsic gluino enhancements at high xp over the leading-twist pre- 
dictions imply that this region of phase space could be especially appropriate for i?-hadron 
searches in the light gluino scenario. For rUg in the 1 — 5 GeV range, a mass region where 
substantial evidence for the analogous intrinsic heavy quark states exists and for which our 
computational techniques should be most reliable, the enhancements are very significant 
(factors of several hundred to several thousand being common). The magnitudes we predict 
for these enhancements may even be conservative since the increased color factor associated 
with intrinsic gluinos compared to intrinsic charm has been neglected. 
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FIGURES 





FIG. 2. The x distribution of intrinsic i?-hadrons in the proton with rUg = 1.2 GeV. The curves 
are g (sohd), W (dashed), 5° (dot-dashed), and (dotted). 
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FIG. 3. Intrinsic gluino higher twist contributions to duig/xp in i?-hadron production with 
i = 1.2 GeV. The solid curve is i?"*", the dotted curve is and the dashed curve is S^. Each 



distribution includes the contribution from independent uncorrelated fragmentation of a gluino. 
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FIG. 4. i?+ xp distribution from 800 GeV protons on various targets. Four gluino masses are 
chosen, rrig = 1.2 GeV (top), 1.5 GeV, 3.5 GeV, and 5.0 GeV (bottom). For each mass, there is a 
triplet of curves representing different targets: proton (sohd), Be (dashed), and Cu (dotted). 
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FIG. 5. xp distribution from 800 GeV protons on various targets. Four gluino masses are 
chosen, rUg = 1.2 GeV (top), 1.5 GeV, 3.5 GeV, and 5.0 GeV (bottom). For each mass, there is a 
triplet of curves representing different targets: proton (sohd), Be (dashed), and Cu (dotted). 
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FIG. 6. R xp distribution from 800 GeV protons on various targets. Four gluino masses are 
chosen, rUg = 1.2 GeV (top), 1.5 GeV, 3.5 GeV, and 5.0 GeV (bottom). For each mass, there is a 
triplet of curves representing different targets: proton (sohd), Be (dashed), and Cu (dotted). 
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FIG. 7. Intrinsic gluino enhancement to xp distribution for various i?-hadrons with rrig = 1.2 
GeV. The lower curve is the fusion basehne for gluino production with a delta function fragmen- 
tation. The curves are (solid), (dashed), (dot-dashed), and the leading twist gluino 
production (dotted). 
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FIG. 8. Intrinsic gluino enhancement to xp distribution for various i?-hadrons with rrig = 1.5 
GeV. The lower curve is the fusion basehne for gluino production with a delta function fragmen- 
tation. The curves are (solid), (dashed), (dot-dashed), and the leading twist gluino 
production (dotted). 
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FIG. 9. Intrinsic gluino enhancement to xp distribution for various i?-hadrons with rrig = 3.5 
GeV. The lower curve is the fusion basehne for gluino production with a delta function fragmen- 
tation. The curves are (solid), (dashed), (dot-dashed), and the leading twist gluino 
production (dotted). 
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FIG. 10. Intrinsic gluino enhancement to xp distribution for various i?-hadrons with nig = 5.0 
GcV. The lower curve is the fusion baschne for gluino production with a delta function fragmen- 
tation. The curves are ii+ (solid), 5° (dashed), (dot-dashed), and the leading twist gluino 
production (dotted). 
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